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We examine the problem of a repulsive, weakly-interacting and harmonically trapped Bose-Einstein 
condensate under rotation. We derive a simple analytic expression for the energy incorporating 
the interactions when the angular momentum per particle is between zero and one and find that 
the interaction energy decreases linearly as a function of the angular momentum in agreement with 
previous numerical and limiting analytical studies. 
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One of the questions that followed the experimental 
realization of Bose-Einstein condensation |j in vapours 
of alkali atoms Q is how they behave under rotation. 
The experimental observation of vortex states in a two- 
component system has been reported by Matthews et al. 
|3j, and the existence of vortex states in a stirred one- 
component Bose-Einstein condensate has been demon- 
strated by Madison et al. Theoretical studies have 
also addressed this problem both in the Thomas-Fermi 
limit of strong interactions between the atoms and 
in the limit of weak interactions [p[-4l3[ . 

In the limit of weak interactions and for a given an- 
gular momentum L, there is a large degeneracy corre- 
sponding to the ways of distributing L units of angular 
momentum among N harmonically confined atoms. Di- 
agonalizing the interaction within this space of degener- 
ate states, Bertsch and Papenbrock (T^] determined the 
ground state energy numerically for a limited number of 
particles and an effective repulsive interaction between 
them; the same method can evidently provide the com- 
plete excitation spectrum of the system. One remarkable 
observation is that, for 2 < L < N, the ground-state en- 
ergy appears to have a particularly simple form with an 
interaction energy which decreases linearly with L. This 
result has also been obtained analytically in Ref. |L2| in 
the region L < N with use of a Bogoliubov transforma- 
tion, and it is also in agreement with the study of this 
problem using the mean- field approximation PJl2]|. 

Here, we present an analytic derivation of this linear 
dependence of the ground state energy on L. Our starting 
point is the hamiltonian H = Hq + V, where 



Ho 



2M 



(1) 



includes the kinetic energy of the particles and the po- 
tential energy due to the trapping potential, and 



S(Ti 



(2) 



is the interaction energy, which we assume to be a contact 
interaction. Here, M is the mass of the atoms, lu is the 



frequency of the trapping potential, which is assumed to 
be isotropic, and Uq — iTrTi 2 a/M is the strength of the ef- 
fective two-body interaction, with a being the scattering 
length for atom-atom collisions. We assume that a > 0, 
i.e., we treat only the case of repulsive interactions be- 
tween the atoms. Further, we restrict our attention to 
the domain 



nUo <C hui, 



(3) 



where n is the density of the condensed atoms. If the sys- 
tem has angular momentum L, this condition allows us 
to work within the subspace of states which are degener- 
ate in the absence of interactions. All other states differ 
by an energy of order hu), which is much larger than uUq 
by assumption. 

It is convenient to work in the basis 



{NoiNuN*,...) = \0 N °,1 N \2 N *,...), 



(4) 



where N k is the number of particles with angular momen- 
tum Tik. These occupation numbers must be chosen to 
respect the obvious constraints on particle number and 
angular momentum: 



(5) 



We introduce annihilation and creation operators and 
at , which destroy or create one particle with angular mo- 
mentum hk. The number operator N and the angular- 
momentum operator L can then be expressed as 



N 



] a\ak] L = 7^ Tika\dk- 
k 



(6) 



We also introduce the raising and lowering operators, L±, 
as 

L+ = Vfc + 1 a\ +l ak] L- = ^ Vk + la\a k+ i. (7) 

k k 

Finally, we write the interaction V, Eq. (||), as 



1 



V 



K 



at a K _ k , 



/ ^ 



Jk 



X 



(8) 



\fc=0 



L k \JlK-k J 

where Ik — kl and Jk — Kl/2 K . It is readily seen that 

{L,L±]=L ± , (9) 

which shows that the operators L± raise or lower the an- 
gular momentum by one unit when applied to the states 
of Eq. (^). Further, both the interaction, V ', and the 
hamiltonian commute with L±, 



[V,L ± ] = [H,L ± }=0. 



(10) 



This shows that, if \L, N) is an eigenstate of V with eigen- 
value £l,n> the state L + \L, N) will also be an eigenstate 
of V with the same eigenenergy but with L + 1 units of 
angular momentum. 

We now proceed to the demonstration that the spec- 
trum is indeed linear if L < N. For any given eigenstate 
with a given N and L, the repeated application of L + 
will generate infinitely many new eigenstates with angu- 
lar momenta L+l, L + 2, .... These states are obtained by 
adding angular momentum to the center of mass coordi- 
nate |p|, ^0| , ^4t while leaving the wavefunction unchanged 
in the relative coordinates of the particles. Since the 
interaction energy depends only on relative coordinates 
and since the center of mass variables separate for a har- 
monic potential, the interaction energy is the same for 
each of these states. For any given N and L, it is useful 
to pay special attention to those "intrinsic" eigenstates, 
\L,N) i n t, which do not have any center of mass excita- 
tion. The observation that the interaction energy in the 
ground state decreases with increasing L [[ll| indicates 
that the ground state wavefunction is an intrinsic state. 
When L + acts on any eigenstate \L— 1, N) (i.e., intrinsic 
or not), it creates a state \L, N) with center of mass exci- 
tation. Such states are necessarily orthogonal to intrinsic 
states with angular momentum L, and thus 

int {L,N\L+\L- 1,N) = (L- l,N\L_\L,N) iat = 0. (11) 

Since the states \L — 1, N) form a complete set with an- 
gular moment L — 1, we conclude from Eq. (11) that 



L_|L,AOint = 0. 



(12) 



Let us now consider two particular basis states, \A) 
and \B), for any L < N in the occupation number repre- 
sentation of Eq. (||) , 

\A) = \N-L + 1,L-2,1); \B) = \N — L, L, 0). (13) 

States | A) and \B) are not eigenstates of V but will in 
general be components of these eigenstates. Similarly, 
consider the state |C) 



\C) = \N-L + 1,L- 1,0). 



(14) 



Note that these three states have = for all k > 3. 
Under the action of L_, \A) and \B) are the only states 
which can lead to state \C). Therefore, if an intrinsic 
eigenstate \L, N) int includes these basis states in the com- 
bination a I A) + (3\B), we see from Eq. (Il2|) that 



L-(a\A) + p\B)) oc \C) = 0. 



(15) 



We can satisfy Eq. (|15|) cither by choosing a = [3 = or 
by insisting that 







2(L-1) 



L(N-L + 1) 



1/2 



(16) 



The second option is interesting since it uniquely deter- 
mines the ratio of two components in \L, iV); nt . 

Since the state \L, N)i nt is an eigenstate of V, it follows 
that 



(B\V\L,N) int = £ LjN (B\L,N) 



ini • 



(17) 



However, the only component of the state \L, N)i nt which 
can be connected to state \B) by the interaction V is 
precisely the linear superposition a\A) + j3\B). Since the 
ratio P/a is known, this fact can be used to determine the 
eigenvalue £l.n- Since \A) and \B) involve particles with 
k = 0, 1, and 2 only, the relevant part of the interaction 
V is thus 

V = l;U (alala ao + 2a\a\a x aQ + ^a\a\aiax 

H — -=a\a\a 2 a n H — —alal^ai). (18) 
V 2 V 2 

The required matrix elements are 

/3 1 

— (B\ ajajoofflo + 2a|ajaia + -a\a\aiai\B) = 

^[(N - L)(N — L — 1) + 2L(N - L) + \h{L - 1)], (19) 



and 
a 
2V2 



B\ ala\a 2 a \A) 



2V2 



y/(N-L-l)(L-l)L. (20) 



It now follows directly from Eqs. ©, @, ©, and © 
that 



N{2N~L~2) 
£l,n — : vq, 



(21) 



in agreement with the numerical results of Refs. [pLl| , |l2[| . 
Here vo = Uo J |$o(r)| 4 dr, and $0 is the ground state 
wavefunction of the harmonic oscillator. To be precise, 
we have demonstrated that an intrinsic eigenstate either 
has the eigenvalue Eq. (|2l] ) or does not contain the com- 
ponents I A) and \B). Completeness indicates that at 



2 



least one intrinsic eigenstate with eigenvalue given by 
Eq. (|2l| ) must exist. 

As we mentioned earlier, the interaction energy of 
Eq. (|2^) varies linearly with L for the contact two-body 
potential of Eq. (Q). This feature is more general. As long 
as one deals with a two-body potential that conserves the 
angular momentum and commutes with L±, this linear- 
ity persists, with the slope 8£l,n / dL being dependent 
on the specific form of the interaction. 

We have exploited the special features of harmonically 
trapped, weakly interacting bosons to determine a sin- 
gle eigenvalue of V from the properties of two (in gen- 
eral exponentially small) components of the wavefunc- 
tion. The present elementary argument does not enable 
us to construct the full eigenfunction or to demonstrate 
that this eigenstate describes the ground state of the sys- 
tem. While the present method can be applied to other 
choices of V, there appears to be no simple generalization 
to the interesting case L > N. 
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